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1. Introduction 




Electroweak gauge boson production with subsequent leptonic 


decays plays a crucial role 


in the physics studies performed at hadron colliders. With larg 


;e cross sections and a clear 



collider signature, the Drell-Yan processes have become a standard candle at the LHC and 
many precise measurements are currently being performed [1-4]. Their study can be used 
to help calibrate detectors and place constraints on parton distribution functions (PDFs), 
see e. g. [5-7]. It has also been proposed that measured electroweak boson production cross 
sections can be used as a luminosity monitor [8]. Drell-Yan processes also constitute an 
important source of background for searches for new physics, such as Z' and W boson 
production and other high mass dilepton resonances [9-12]. 

Given the large amount of data that the LHC will deliver, measurement errors will 
eventually be dominated by systematics rather than statistics, see e. g. [13]. Consequently, 
electroweak boson production needs to be known to high precision. Much effort has been 
made in this area. Indeed, it was the first hadronic scattering process to be computed at 
next-to-next-to-leading (NNLO) in QCD [14,15]. Now also the differential cross sections 
for Drell-Yan processes are known through NNLO in QCD [16-21], and the NLO QCD cor- 
rections have been matched to parton showers [22-25]. Electroweak corrections have been 
known for sometime [26-29], and continue to be refined [24,30-34]. Attempts to approxi- 
mate QCD X EW corrections to Drell-Yan have also appeared in the literature [34,35] and 
very recently the mixed QCD x EW two-loop virtual corrections have been calculated [36]. 



- 1 - 



At hadron colliders, the high center-of-mass energy (cm.) typically leads to produc- 
tion of electroweak gauge bosons in association with QCD radiation. Here, we focus on 
Drell-Yan processes with one additional, hard jet. The intermediate gauge boson then 
recoils against the jet and can be strongly boosted. If the jet has large transverse momen- 
tum, the leptons originating from the gauge boson decay are produced at high transverse 
momentum as well and provide a source for either high-energy dilepton pairs of opposite 
charge (Z/7*+jet production) or high missing transverse energy in combination with a 
charged lepton (Pl^+jet production). Theoretical predictions for the Drell-Yan+jet cross 
section include NLO QCD corrections [37-39] and parton shower matching [22], as well as 
in the electroweak sector the NLO corrections in the on-shell approximation [40-44] and 
the full NLO calculation, including photonic corrections and including the leptonic decay 
of the gauge boson [45,46]. Also for Drell-Yan cross sections with higher jet multiplicities, 
the NLO QCD corrections are now available [39,47-50]. 

Often, precision measurements of Standard Model (SM) processes are used to constrain 
new physics, not by direct detection, but by their influence through quantum fluctuations. 
One may ask the question how radiative corrections from physics beyond the Standard 
Model (BSM) would effect such well-known processes. To allow for precision tests of the 
SM and its perturbative expansions, it is crucial to know if there were large BSM corrections 
since otherwise information about the underlying physics cannot reliably be extracted from 
the data. An interesting paradigm for physics beyond the SM is supersymmetry (SUSY). 
Many collider studies have been performed to examine the phenomenology of SUSY pro- 
duction at hadron colliders, and the possibility of their detection. The investigation of 
new physics contributions to SM processes via higher-order corrections is a logical next 
step. Only in recent years, Drell-Yan processes have been studied within the framework 
of the minimal supersymmetric Standard Model (MSSM). The NLO supersymmetric QCD 
(SQCD) and electroweak corrections within the MSSM to the charged- and neutral-current 
Drell-Yan processes with no final state parton were calculated in [30,31] and the NLO elec- 
troweak corrections for the on-shell VF+jet production process within MSSM have been 
computed in [51]. The impact of SUSY corrections has been found to be small. 

In this work, we further complete the one-loop picture and calculate the NLO SQCD 
corrections within the MSSM to Drell-Yan processes in association with a hard jet in the 
final state. We consider both Z'/7*+jet and VF+jet production at the LHC, including the 
leptonic decay of the gauge bosons, and take all off-shell effects due to the finite widths of 
the Z and W boson and all contributions of an intermediate photon into account. Our aim 
is to show the stability of Drell-Yan+jet processes under the inclusion of SQCD corrections 
within the MSSM. 

This paper is organized as follows. The setup of the calculation is detailed in section 2. 
The computation of the SQCD corrections to the neutral- and charged-current Drell-Yan 
processes proceed in close analogy and are discussed in parallel. In section 3 we present the 
numerical results for both Z/7*+jet and VF+jet production at the LHC. We summarize our 
findings in section 4. Details on the renormalization procedure are given in appendix A. 
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Figure 1: Feynman diagrams for the quark-anti-quark initial state process, Eq. (2.1). 
2. Details of the calculation 

In this section we describe the technical setup of our calculation of SQCD corrections for 
the neutral- and charged-current Drell-Yan process with an associated hard jet. We always 
include the decay of the gauge bosons into leptons and take all off-shell effects due to the 
finite width of the Z ov W boson into account. We first address the leading order (LO) 
processes in section 2.1 and then discuss the NLO calculation in section 2.2. In view of 
SQCD corrections, the computations for Z/j*+iet and VF-l-jet production are very similar 
and proceed in close analogy. 

2.1 Leading order processes 

At the LHC, the Drell-Yan process in association with a hard jet is initiated by three differ- 
ent partonic processes at LO. Including the leptonic decay of the intermediate electroweak 
gauge boson, they are in the case of the neutral-current Drell-Yan process, 

qq^ Z/^* g^re-^ g , (2.1) 
g q ^ q ^ ri+ q , (2.2) 

gq^ Zif q ^ r ^+ q , (2.3) 

and for the charged-current Drell-Yan process, 

q(i g^lv.g , (2.4) 

gq^W (i ^tvicl , (2.5) 
gq^W 4 ^ ■ (2.6) 

For production, a positively charged lepton is produced together with an anti-neutrino 
{t^vi)^ for production it is a negatively charged lepton and a neutrino {t~V(). 

The tree-level Feynman diagrams for the process Eq. (2.1) can be found in Figure 1, 
those for l^-|-jet production follow analogously by replacing the Zqq (Zii) vertex with 
a Wqq' iWlvf) vertex. Diagrams for the (anti-)quark-gluon processes are obtained by 
crossing the gluon with an initial-state (anti-) quark. 

The quarks considered are q, q' = u, c, d, s with q and q' being quarks of opposite 
isospin. The quarks are treated as massless throughout the calculation. The Z/7*-|-jet LO 
cross section is quark generation independent, aside from the appropriate PDF inclusion 
and convolution. The only up- and down-type quark dependence comes from the vector 
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and axial couplings in the Zqq vertex. The VF+jet cross sections do depend on the quark 
flavor due to the quark mixing parameterized by the CKM matrix. However, the CKM 
quark mixing factorizes from the tree-level matrix elements and, for the inclusive cross 
section, there needs only one squared amplitude to be calculated for each of three partonic 
channels, weighted by the sum of squares of the respective absolute value elements of the 
CKM matrix, as well as the corresponding PDFs. 

Throughout the calculation, we exclude bottom-quark initiated processes due to their 
PDF-suppressed small contribution to the cross section and we neglect the mixing between 
first and third generation quarks. A fixed width is included in the propagator of the 
intermediate Z or W boson, giving the expected Breit-Wigner resonance. 

Already at LO, there exists the possibility that the parton in the final state may become 
soft or collinear to the gluon or quark that it couples to. If this occurs, the cross section 
becomes infinite. Divergences of this type are regulated by requiring the jet to be hard. 
This is achieved by applying a minimum transverse-momentum (pt) cut on the jet. Also 
we treat the leptons as massless throughout the calculation and require a minimum lepton 
Pt as well as a minimum lepton-pair invariant mass, M^, in case of Z/7*-|-jet production. 

2.2 NLO SQCD corrections 

We calculate the one-loop SQCD corrections to the partonic processes in Eqs. 2.1-2.3 and 
Eqs. 2.4-2.6, which are purely virtual. The SUSY particles charged under SU{3)c are 
the squarks, iji, and gluinos, g, the superpartners of the SM chiral quarks and gluons, 
respectively. In the SM, left- and right-handed quarks transform differently under SU{2)l- 
Consequently, under supersymmetric transformations, there are separate superpartners, 
Ql, Qb.i for quarks of different handedness, ql, qr. After electroweak symmetry breaking 
(EWSB), left- and right-handed squark eigenstates mix to form squark mass eigenstates, 
qi, q2- Squark mixing, however, is proportional to the mass of their SM quark partner. 
In this calculation, only squarks of the same fiavor as the initial state quarks are present. 
Since the masses of the quarks considered (u, d, c, s) are small, we can neglect any left-right 
mixing between squark eigenstates, and the two gauge eigenstates of a given flavor q are 
also mass eigenstates that we denote hy iji, i = L, R. 

In the calculation of Z/7*-|-jet production, the relevant interactions between SM and 
SUSY particles are the 'yqiijj, Zqiq.j, gtjiqi, ggg, and qtjig vertices. Although there are 
diagrams with the ^gqiqj or Zgqiqj vertex, their contribution to the virtual corrections is 
zero, and they have not been explicitly included in this manuscript. There are 44 diagrams 
per partonic channel that give non-zero contributions to the differential cross section (22 
for Z boson and photon mediation each). They divide into self-energy insertions, gluon- 
and Z/7*-vertex corrections, and box contributions and are organized in Figure 2. In 
each diagram, the two squark eigenstates qi of the same quark flavor q as the initial state 
can run in the loop. The Feynman diagrams for l/F-|-jet production can easily be inferred 
from Figure 2 by again replacing the Zqq (Zii) vertex with a Wqq' (Wii'e) vertex and 
by replacing the Zijiiji by Wqiq[ vertices. Of course the W boson couples to left-handed 
particles only, and only the left-handed components of the squark eigenstates contribute 
(i. e. qi = qi in case of no left-right mixing). 
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(a) Self-energy insertions: 




(b) Gluon-vertex corrections: 




(c) Z/7*-vertex corrections: 




(d) Box contributions: 




Figure 2: Feynman diagrams for tiie self-energy, gluon- and Z/7*-vertex, and box contributions 
to the partonic process qq i~^£~g, mediated by Zf-y* excliange, at SQCD NLO. 

Due to only massive particles propagating in the loop, the virtual corrections are com- 
pletely infrared (IR) safe. Therefore, there are no real emission processes to be considered 
to cancel singularities. However, as noted above, we require a hard jet in the final state. 
This statement implies a minimum cut has been placed on the transverse momentum of 
the jet to avoid IR divergences that would otherwise appear at LO. 

Ultraviolet (UV) divergences arise from the self-energy insertions and the gluon and 
Z/'j* or W vertex corrections, while the box diagrams are IR and UV finite. Dimensional 
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regularization is used to regulate the UV divergences, where the number of dimensions is 
d = 4 — 2e and the singularities then take the form of 1 /e poles. In order to remove the UV 
divergences one has to include the proper counterterms for SQCD one-loop renormalization. 
Our renormalization procedure closely follows ref. [52]. We impose on-shell conditions to 
fix the renormalization constants in the quark sector and use the MS scheme, modified to 
decouple the heavy SUSY particles [53,54], for the renormalization of the strong coupling 
constant and the gluon field. The explicit expressions in terms of the renormalization 
constants can be found in appendix A. 

We set up the calculation in the conventional diagrammatic approach. As an im- 
portant cross-check on the numerical results, two completely independent computations 
were performed using different techniques. The first calculation is based on the framework 
of FeynArts 3.6 [55], FormCalc 7.0, and Loop Tools 2.6 [56]. In the second calculation, 
Feynman diagrams are generated with QGRAF [57]. Algebraic manipulation is performed 
with an in-house software comprised of Maple and Form [58]. Then AIR [59] is used for 
the reduction of tensor loop integrals to scalar ones. The numerical integration is per- 
formed using Vegas of the Cuba 1.7 [60] library, as well as the scalar loop integral library 
QCDLoop [61]. 

The Zqq, Wqq', Zqiijj, Wqiq'p and qqg vertices contain axial pieces. Special attention 
should be paid to the the treatment of 7^ when working in d dimensions. In the first 
calculation the naive anticommuting scheme [56, 62] was used. Here, it is assumed that 
{7^)7^} = as in four space-time dimensions and that the four-dimensional, non-zero 
trace, 

Tr[7^7^7V7^] = -4ie^^'"", (2.7) 

remains. This FormCalc-based approach has been justified in [63]. In the second calcula- 
tion, the treatment of 7^ as described in [64] is used, where 

7' = ^e^.p.l^YYr ■ (2.8) 
Both approaches have been shown to give consistent results in our NLO calculation. 



3. Numerical results 

In this section we investigate the numerical influence of the SQCD one-loop corrections to 
both the neutral- and charged current Drell-Yan process with an associated hard jet. We 
examine these processes at the LHC with a c. m. energy of 7 TeV, unless noted otherwise. 
We list the relevant input parameters in section 3.1. The integrated cross section results 
and differential distributions for the Z/7*-|-jet process are given in sections 3.2 and 3.3, 
respectively. The results for l^-l-jet production are discussed in sections 3.4 and 3.5. 
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3.1 Input parameters 

The SM input parameters are chosen in accordance with [65], 

Mz = 91.1876 GeV, Mw = 80.339 GeV, 

Tz = 2.4952 GeV, Tw = 2.085 GeV, 

cos^w^ = Mw/Mz, a-^ = aiMz)"^ = 128.91, ^^'^^ 

iKdl = 0.97428, |K.| = 0.2253, l^dl = 0.2252, |Fc.| = 0.97345. 

Quarks and leptons are considered massless and we do not further specify the flavor of 

the leptons in the final state as the SQCD corrections do not depend on the lepton flavor. 
Factorization and renormalization scales are fixed and identified with the vector boson 
mass, ij,F = m = Mz,w- We use the central MSTW2008NLO (68% CL) PDF set [66] in 
its LHAPDF implementation [67], with the strong coupling OsifiR) they provide, yielding 

asiMz) = 0.12018, asiMw) = 0.12257. (3.2) 

We require a set of basic kinematic cuts to be satisfied. As previously mentioned, a 
minimum pT of the final state parton is required to render the cross section finite. We also 
demand the two final-state leptons to have a minimum transverse momentum. Furthermore 
we require the leptons and the jet to be produced in the central region of the detector and 
apply a cut on their rapidities. For the Z/^*+jet process we also require a minimum 
lepton-pair invariant mass, M(£. We choose the following numerical values. 



Z/7*+jet : pr,jet > 25 GeV, |yjet| < 2.5, 

PT,e± > 25 GeV, |y^± | < 2.5, Mu > 50 GeV . 



(3.3) 



t^+jet : pT,iet > 25 GeV, |yjet| < 2.5, 
Pr,^ :?^r > 25GeV, |y^| < 2.5 . 

However, only the integrated cross sections depend strongly on the specific cuts chosen, 
while we have found that the relative importance of the SQCD corrections does not vary 
much when the cuts arc tightened or loosened. 

The only SUSY particles that appear in the loops are squarks and the gluino. The 
SQCD corrections are fiavor- and chirality-blind and no other SUSY parameters than the 
squark and gluino masses enter the calculation. For simplicity we neglect the squark left- 
right mixing, set all squark masses equal and use a common squark mass, rriq, and the 
gluino mass, rrig, as direct input. There is no need to define a complete set of SUSY input 
parameters and we do not consider commonly used benchmark scenarios (as e. g. SPSla') 
here. As we will see below, this approach is sufficient for the purpose of our study to 
show the stability of Drell-Yan-|-jet under the inclusion of SQCD corrections. We use the 
following values for our numerical studies, if not otherwise noted, 

my- = 600 GeV, rrig = 500 GeV. (3.4) 

These sparticle masses are already at the lower limit of the mass region that is currently 
investigated by LHC experiments, see e. g. [68,69], and allow for a conservative estimate of 
the typical size of SQCD corrections to Drell-Yan-|-jet processes. 
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Z/7*+jet 


partonic 


LO 


SQCD 


5 


production 


channel 


cross section 


contributions 


ys = / lev 




16.55 pb 


xj.yoi ID 


U.UUoy /o 




fin -\- on 


37 78 nb 


2 52(1 fb 


0066 % 




incl. 


54.27 pb 


3.507 fb 


0.0065 % 


= 14TeV 




31.70 pb 


1.947 fb 


0.0061 % 




9Q + gq 


91.38 pb 


6.798 fb 


0.0074 % 




incl. 


123.1 pb 


8.746 fb 


0.0071 % 



Table 1: Numerical results for integrated cross sections for the neutral-current Drell-Yan process 
mediated by a Z boson or virtual photon 7* in association with a hard jet at the LHC, with 
y/s — 7TeV and ^/s — 14TeV. Shown are the leading order results in picobarn (pb), the SQCD 
contributions in femtobarn (fb) and the relative corrections 5 for the partonic subchannels and the 
inclusive result (inch). Light quarks are implicitly summed over in the initial state, q — u,d,c,s. 
We consider = 600 GeV, rrig = 500 GcV and the cuts listed in Eq. (3.3) have been applied. 
Factorization and rcnormalization scale are set to /i = Mz (with MSTW 2008 NLO) . 

3.2 Z/7*+jet integrated cross section results 

Here we present the integrated cross section for charged dilepton production with a hard 
jet at the LHC. Table 1 shows the LO cross sections and SQCD contributions for the qq 
and gluon-initiated partonic processes, at a proton-proton c. m. energy of y/s = 7TeV and 
= 14TeV. 

The prominent production modes at the LHC are the gluon-induced initial states, 
enhanced by the large gluon densities at small parton momentum fractions. With the 
kinematic constraints in Eq. (3.3), the total integrated cross section is 54.27 pb at 7 TeV 
and 123.1 pb at 14 TeV {q = u,d,c,s). The SQCD corrections, for squark and gluino 
masses of 600 and 500 GeV, respectively, account for an increase of 0.006% — 0.007% in 
the inclusive integrated cross section as well as per partonic channel. The c. m. energy at 
the LHC has little effect on the size of the corrections. Only the LO cross section strongly 
depends on the energy. When going from 7 TeV to 14 TeV, the gluon-quark channels 
become even more important and we find that at LO they contribute 74%. However, the 
SQCD effects on the integrated cross section remain basically unchanged and small in 
absolute value. 

To further study the effects of SQCD corrections on the integrated cross section of 
Z'/7*-|-jet production, we show in Figure 3 the relative corrections as a function of the 
common squark mass niq, for different values of the gluino mass. As one would expect, the 
impact of the SQCD corrections increases as the SUSY particle masses decrease. However, 
even for very light squarks and gluinos of only 100 GeV the relative corrections are still 
significantly below the 1% level. Figure 3 shows that the Z/7*-|-jet integrated cross section 
is stable under SQCD corrections, over a broad range of low-mass squarks and gluinos. 
For TeV-range SUSY particles the SQCD corrections are completely negligible. This is 



-8- 



0.25- 



Z/y'+jet — 








100 GeV 
200 GeV 
300 GeV 
400 GeV 
500 GeV 



100 



200 



300 
m(q) [GeV] 



400 500 600 



1000 



Figure 3: Relative SQCD corrections of the integrated cross section, S = [ajsiLO ^ '^lo)I'^lOi for 
charged lepton-pair production mediated by a Z boson or virtual photon in association with a hard 
jet at the LHC, with -y/s = 7 TeV, as a function of a common squark mass, m^, for different gluino 
masses, rrig. The cuts listed in Eq. (3.3) have been applied. 

comforting in the sense that if Z/7*+jet is used as a normalization process, no SUSY BSM 
physics should be masked in this normalization. In terms of our calculation, this stability 
also justifies our neglecting of squark left-right mixing, and use of a degenerate squark mass 
scheme. 

3.3 Z/7*-|-jet kinematic distributions 

We have seen above that the SQCD contributions to the integrated cross section are sub- 
tle. However, the corrections can be more pronounced in the differential distributions of 
kinematic variables. We consider the LO and NLO differential cross sections, and define 
the relative corrections 5,5= '^'^^§~q^'^ j foi' given observable O, where Onlo is the 
sum of the LO contributions and the SQCD contributions. We present distributions in the 
lepton-pair invariant mass, M^, the lepton transverse momentum and rapidity, and 
yi, and the jet transverse momentum and rapidity, PTjet and yjet- All distributions shown 
have been subject to the kinematic constraints given in Eq. (3.3). 

In Figure 4, distributions in the lepton-pair invariant mass Mu are displayed. We find 
that the SQCD corrections hardly affect the shape of the LO distribution and are below the 
1% level for Mu < 1 TeV (left panels). In the vicinity of the Z boson resonance, shown in 
the upper right panel, they are completely negligible and do not distort the SM result. For 
larger values of Mu (lower right panel) the relative corrections can reach several percent, 
until a SUSY mass threshold is reached and then the relative corrections begin to fall. 

The lepton differential distributions in transverse momentum, pt/, and rapidity, y£, 
are given in Figure 5. In the low- and intermediate-py range, the SQCD contributions to 
the pT,i distribution are positive but small, increasing only to a maximum of about 1% 
at around the threshold of the average sparticle mass. In the high-p^^ region, the relative 
corrections become non- negligible, on the order of a few negative percent. Concerning the 



-9- 




Figure 4: Differential distributions for the neutral-current Drell~Yan process in assocation with a 
hard jet at the LHC. NLO SQCD corrections are calculated for niq = 600 GeV and nig — 500 GeV 
and the cuts given in Eq. (3.3) have been placed. Shown are the absolute differential distributions for 
the LO and NLO processes (top left) and the relative difference between NLO and LO distributions 
(bottom left) with respect to Icpton-pair invariant mass Mu. In the right panels, the relative 
corrections around the Z boson resonance and in the higli-M^^ region are shown, respectively. 



lepton rapidity. Figure 5 (right), we find that the relative corrections are nonzero but far 
below 1% and thus do not have any impact on the shape of the SM-only distribution. It 
is important to observe the discrepancies between the lepton and anti-lepton distributions 
in Figure 5. These differences occur already at the LO process and are a result of the Z 
boson's axial coupling to fermions. Setting the axial coupling to zero would render the 
lepton and antilepton curves identical. 

The jet distributions. Figure 6, exhibit the same behavior as the lepton distributions. 
We see the threshold effect around {niq + mg)/2 is even more pronounced and the relative 
corrections in the high-p^jet region [pT > 1.5 TeV) contribute with roughly negative 4 — 5%, 
Figure 6 (left), while there are vanishing corrections to the jet rapidity, Figure 6 (right). 

We also investigated differential distributions for lighter SUSY particles. We found a 
very similar picture as described above, with the SUSY threshold shifted accordingly, and 
abstain from showing the plots here. For rriq = 400 GeV and rrig = 350 GeV, the SQCD 
corrections become only slightly more important, and reach the 1 — 2% level at the sparticle 
thresholds and the —5% level for px ~ 1500 GeV. 

3.4 VF+jet integrated cross section results 

In Table 2 we show the results for the integrated cross section for lepton-neutrino production 
with a hard jet, mediated by a boson, at the LHC. The LO cross sections and SQCD 
contributions for the qq' and gluon-initiated partonic processes are given, at a proton-proton 
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Figure 5: Lepton transverse momentum, pt,i, (a) and lepton rapidity, y^, (b) distributions for the 
Z/7*+jet process at the LHC with ^/s — 7 TeV. The black line represents the lepton while the 
red line represents the antilepton. NLO SQCD corrections are calculated for rriq — 600 GeV and 
nig = 500 GeV and the cuts given in Eq. (3.3) have been placed. 




Figure 6: Jet transverse momentum, PTjetj (a) and jet rapidity, ?/jct, (b) distributions for the 
Z/7*+jet process at the LHC with ^/s = 7 TeV. NLO SQCD corrections are calculated for nig = 
600 GeV and nig — 500 GeV and the cuts given in Eq. (3.3) have been placed. 
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48.52 fb 
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1056 pb 


59.33 fb 
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Table 2: Numerical results for integrated cross sections for the charged-current Drell-Yan process 
mediated by a boson in association with a hard jet at the LHC, with — 7TeV and y^s = 
14TeV. Shown are the leading order results in picobarn (pb), the SQCD contributions in femtobarn 
(fb) and the relative corrections 5 for the partonic subchannels and the inclusive result (incl.). Light 
quarks are implicitly summed over in the initial state, q = u,d, c, s. We consider rriq = 600 GeV, 
rrig — 500 GeV and the cuts listed in Eq. (3.3) have been applied. Factorization and renormalization 
scale arc set to fi ^ Mw (with MSTW 2008 NLO) . 

cm. energy of ^/s = 7TeV and ^/s = 14TeV. Similar to the Z/7*+jet case, the SQCD 
corrections contribute between 0.004 — 0.006% to the individual channels and the inclusive 
result. The gluon- induced processes are slightly more affected than the quark-annihilation 
channels. The main cross section contribution comes from gluon-induced initial states, with 
increasing importance for higher c. m. energies. Thus the impact of the SQCD corrections 
on the inclusive cross section is somewhat enhanced at 14TeV compared to 7TeV. 

3.5 Vl^+jet kinematic distributions 

The differential distributions for VK+jet production at the LHC with ^/s = 7TeV are shown 
in Figs. 7, 8, and 9. They correspond to distributions in transverse mass, lepton transverse 
momentum and rapidity, and jet transverse momentum and rapidity. 

The transverse mass, Mt{(.v)^ and the pxi distributions are particularly relevant for 
the measurement of the W boson mass at hadron colliders. Here, the transverse mass is 
defined as Mt{Iu) = [{\ptA + \i>T\? " (Pt,^ + M^Y'""- 

Again, the results are very similar to those for Z/7*+jet production. The MTi^y) 
distributions, displayed in Figure 7, receive percent-level corrections due to SQCD effects, 
which are minimal and almost vanishing when the W boson is on-shell and maximal around 
squark and gluino thresholds in the high- My (^z/) region. 

The relative corrections to the pt/ and PTjct distributions, Figure 8 (left) and Figure 9 
(left), peak around the average sparticle mass where they amount to about 1%, and grow 
negative in the high-pr region, reaching a couple percent in the TeV region. The SQCD 
corrections in the lepton and jet rapidity distributions, Figure 8 (right) and Figure 9 (right), 
are flat and can safely be neglected in an experimental analysis. 
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Figure 7: Differential distributions for the charged-current Drell-Yan process in assocation with a 
hard jet at the LHC. NLO SQCD corrections are calculated for niq = 600 GeV and nig — 500 GeV 
and the cuts given in Eq. (3.3) have been placed. Shown are the absolute differential distributions for 
the LO and NLO processes (top left) and the relative difference between NLO and LO distributions 
(bottom left) with respect to transverse mass Mx{£i')- In the right panels, the relative corrections 
around the W boson resonance and in the high- MT^iiy) region are shown, respectively. 




Figure 8: Lepton transverse momentum, pt,i, (a) and lepton rapidity, yi, (b) distributions for the 
M^+jet process at the LHC with -y/s = 7 TeV. The red line represents the charged lepton while the 
black line represents the neutrino (missing transverse momentum, ^y). 



- 13 - 




200 400 600 800 1000 -3 -2 -1 1 2 3 

p^(jet) [GeV] y (jet) 



Figure 9: Jet transverse momentum, PTjct, (a) and jet rapidity, ?/jct, (b) distributions for the 
W+iet process at the LHC with y/s = 7 TeV. 

4. Conclusions 

In this manuscript, we present a full treatment of the SQCD corrections to the neutral- and 
charged-current Drell~Yan process in association with a hard jet. We include the decay of 
the electroweak gauge bosons into a lepton pair or a lepton-neutrino pair, respectively, and 
take all off-shell effects due to the finite Z and W boson widths into account, as well as 
the contributions and interference effects of the photon-mediated diagrams in case of the 
neutral Drell-Yan process. 

Even though the Drell-Yan process, with an associated jet, can not be used in any 
practical manner to set indirect limits on SUSY particle masses or SUSY parameters, it 
is important to study the impact of SUSY (and other BSM physics) contributions to SM 
processes. As standard candles, and important backgrounds to new physics at the LHC, it 
is necessary to understand the stability of EW gauge bosons production processes against 
BSM contributions. 

We find that the relative corrections to the integrated cross section due to SQCD 
corrections are small, below the 0.5% level even for very light SUSY particles. Examining 
the differential distributions in conventional kinematic variables such as pT, Mu, or MT{ii^), 
we find vanishing effects in the vicinity of the Z or W boson resonance peak, while the 
relative SQCD corrections can increase to 2 — 5% when the electroweak gauge bosons are 
far off-shell. 
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A. Counterterms and renormalization constants 



Here we list the counterterms for the renormalization of vertices and propagators in the 
SQCD one-loop amplitudes for Drell-Yan-|-jet production. All quarks are considered mass- 
less. 

By using multiplicative renormalization, we replace in the QCD Lagrangian the left- 
and right-handed quark fields, ^'l,/?, the gluon field, G", and the strong coupling constant, 
9s, with 

*L,i?^ G^^V^G^, gs^Zg^g,. (A.l) 

Expanding Zi = 1 + SZi, and introducing SZv,a = ^{SZl ± SZr), we find the follow- 
ing Feynman rules for the self-energy and vertex counterterms that are relevant in our 
calculation. 



i^(5Zy - dZA'Jd), 



(A.2) 




: -iesgj'' {5Zv - 5Za1^), 



(A.3) 




cwsw 



r{g'y-g\l^){5Zv-5ZAi^), 



9v ~ 2^3 ^q^W 
9 A — 2-'3' 



' (A.4) 




2\/2sw 



V^I^^^(1-^^){5Zv + 5Za), 



(A.5) 




-igsT^ 1^ {bZy - 6Zai^ + 5Zg^ + Uzg) 



(A.6) 



where Cq is the electric charge of quark q and /g denotes the eigenvalue of the third com- 
ponent of its weak isospin. T", a = 1..8, are the color matrices of SU{?>)c, while we have 
omitted the color indices of the quarks. We use the abbreviation sw and cw for the sine 
and cosine of the electroweak mixing angle 9w- As usual, q and q' are quarks of opposite 
isospin and V^^^^ is the corresponding entry of the CKM quark mixing matrix. 



For the renormalization of the gluon field and the strong coupling we use the MS 
scheme, modified to decouple the squarks and the gluino [53,54]. The renormalization 
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constant of the strong coupling is then given by, see also [52], 

where the sum runs over the twelve squark eigenstates, fi is the renormalization scale, 
A = l/e — 7E + ln47r; and 

/3o = -'^N-^ng, (A.8) 

with = 3 and riq = 6 for the six squark flavors, includes only the gluino and squark 
contributions. In this case, there is a simple relation between the gluon field and strong 
coupling renormalization constants, 

6ZG = -2 6Zg^, (A.9) 

and they actually cancel out in Eq. (A. 6) and do not enter our calculation. 

In the quark sector, we use the on-shell scheme to fix the renormalization constants 
SZy^A of the (massless) quarks, i. e. the renormalization constants are obtained by evalu- 
ating the vector and axial components of the quark self-energy at the on-shell quark mass, 

6Zv,A = -Sy,A(/ = 0). (A. 10) 

The SQCD corrections to the quark self-energy consist of a squark-gluino bubble insertion 
in the quark line and can be expressed in the following general, compact form [52], 

' 2 ^^-^^^ 



i 

Here, , m^) is the two-point function as defined in [56]. i sums over the two 

squark mass eigenstates and and gp are scalar and pseudoscalar couplings. In the limit 

L Ft L R 

of no left-right mixing in the squark sector, g^' = ±1 and gp = 1. When considering 
the mixing, then 

gs' =cos9ggg' ±sm9ggg' , 

1,2 a L,R,.o R,L (.A.iiJ 

g^ =cosdqgp' ±sindqgp' , 

where Og is the squark mixing angle (notation as in [52]). 

In our numerical study, we neglect the left-right squark mixing, and consider degen- 
erate squark masses, = rriq^ = rriq. In this limit, we find S^(p^) = and the only 
renormalization constant that remains is 5Zv, with 



a 



3mi — mi mi I mj,\ m?(2m? — m?) 

— _ A H - - In I - I + — In 

^ 37r 1 2(m|-m|) (m| - 7n|)2 \tJ-pJ ("^| - m|)2 




(A.13) 
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